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Abstract 

< 

f — . We discuss algebraic and analytic structure of rational Lax operators. With algebraic 

reductions of Lax equations we associate a reduction group - a group of twisted automor- 
phisms of the corresponding infinite dimensional Lie algebra. We present a complete study of 
dihedral reductions for si (2, C) Lax operators with simple poles and corresponding integrable 
equations. In the last section we give three examples of dihedral reductions for sl(N, C) Lax 
operators. 
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^ ! 1 Introduction 
O 

The inverse spectral transform method (1ST) for integrating nonlinear partial differential equa- 
tions (PDEs) is based on the Lax representation (see for example pQ, [2], [3]). In this setting, 
the nonlinear equation is equivalent to the compatibility condition of a pair of linear equations, 
called Lax pair. In general, given an integrable nonlinear equation, there is no algorithmic way 
to find its Lax representation. Some progress in the problem of finding a Lax representation for 
a given nonlinear equation have been made by the Wahlquist and Estabrook method of pseudo- 
potentials and prolongation structures [I] , [2] , [S] • In some cases the Painleve approach provides 
both integrability conditions and Lax representation (see for example [7j). 
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A different approach to this problem consists in starting from a quite general Lax repre- 
sentation which yields to a rather big (but integrable) system of equations. In general, such 
a system has too many degrees of freedom and is very complicated. However, it can contain 
smaller subsystems, relevant for the applications: many famous integrable equations (such as 
the sine-Gordon equation and the Tzetseica equation) are indeed the result of reductions of 
more general integrable systems. This observation suggests to study and classify possible Lax 
representations and their reductions. 

The first attempt to study algebraic reductions of Lax representations has been made in [H] , 
[S], and later in The principal observation in the papers [H], |5| was that reductions 
of integrable equations can be associated with a discrete symmetry group (called reduction 
group) of the corresponding linear problems. The simplest example of such a symmetry is the 
conjugation for self-adjoint operators. 
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In this paper, which is meant to be the first of a series, we report new developments of 
this approach. The paper is organised as follows. In Section [5] we revisit the concept of Lax 
representation and its reductions. We define a quite general family of Lax operators with 
rational dependence on the spectral parameter A, on a simple Lie algebra q. In Section |3] wc 
study reductions of integrable equations corresponding to rational Lax pairs, with simple poles 
in I and g = s/(2,C). Section concerns with the study of possible reductions a fairly general 
Lax pair whose underlying algebras is q = sl(N,C). In both sections, reductions are obtained 
using realisations of the reduction group isomorphic to the dihedral group Bat. Some of the 
systems of equations obtained in these sections are new, to the best of our knowledge. Section 
is devoted to conclusions and remarks. 



2 Lax representation of integrable systems 

Let us start this preliminary section reviewing the concept of Lax representation of an integrable 
system. Here and afterwards subscripts mean partial differentiation. 



2.1 Example of a rational matrix Lax pair and gauge transformations 

Let us consider two linear equations 

I4 = ij> x -Q(x,t)il> = 0, M^j = V>t - P(x,t)i> = 0, (1) 

where tp, Q and P are N x N matrix functions of x and t. This is an over-determined system 
of equations and its compatibility condition 

Qt-P x + [Q,P] = (2) 

is equivalent to a system of nonlinear partial differential equations (PDEs). Its general solution 
is given by 

Q = 9x9~ x , P = gtg' 1 , 

where g = g(x, t) is any nonsingular matrix function of x and t. The same problem becomes 
nontrivial if one assumes that the matrices Q and P depend also on an auxiliary parameter 
A, called spectral parameter, and are polynomial or rational functions of this parameter. As a 
consequence, equation @ has to be satisfied for all values of A and this requirement leads to a 
nontrivial system of integrable, partial differential equations. In the latter case, one has to solve 
a certain matrix Riemann-Hilbert problem to find a suitable matrix function g(x,t,X), which 
would provide the solution of the system (see for example pQ, ^2], JH]). In order to receive a 
nontrivial set of nonlinear equations integrable by the 1ST method an essential dependence of 
Q and P in @ on the spectral parameter A is required, where essential dependence on A means 
that it cannot be eliminated by a A dependent gauge transformation (see below). 
As an example, let us consider the case of Q and P given by 

n=l n=l 

where it is assumed that the sets of complex constants {a±, ...,ajVi|On 6 C} and {b%, ...,b]y 2 \b n £ 
C} have empty intersection and the matrix-coefficients Q n ,Pn may depend on x and t only. 
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Here and hereafter I 6 C. Under these assumptions the compatibility condition @ leads to a 
system of N\ + iV2 + 1 equations 



Qot — Pox + [Qo, Po] — 

^ r r-. X ^ ? 



(4) 



(5) 




(6) 



for JVi + N2 + 2 functions, i.e. the system obtained is under-determined. 

Lax representations of this type and their reductions will be studied in Section G3 

The system of equation Q)-© is invariant with respect to transformations 



Q -> g Qg- g g. 



X 5 



P -> g-'Pg - g-'gt 



(7) 



where g is any nonsingular matrix function of x and t (in principle it may depend on the spectral 
parameter A, but in this section we shall assume g = g{x,t)). These transformations are called 
gauge transformations ^2j, In terms of matrices Q n and P n the gauge transformation 
takes the form 

Qo -> g~ l Qog - g~ l g x , Po -> g~ l Pog - g~ l gt , Qn -> g~ l Q n g , P n g~ l P n g , (8) 

where n > 1. There are several ways to fix the gauge. The natural or canonical choice 12 is to 
set Qo = Pq = 0. This is possible because the equations 



are compatible and their compatibility condition coincides with (JIJ. In this gauge, the system 
becomes a well determined system of N\ + N2 matrix equations lfH|) for A^i + N2 matrix 
functions Q n ,P n , n > 1 . Another useful way to fix the gauge is to assume that the matrix 
Qi is in a canonical (for instance diagonal) form. If the eigenvalues of Q\ are distinct, then 
the remaining gauge freedom, i.e. transformations (jSJ) which do not change Qi, consists of 
nonsingular diagonal matrices g. This remaining freedom can be used to make the matrix Qq 
off-diagonal. Such choice of the gauge is called the pole gauge ^2j and it also provides a well 
determined system of equations. 

The nature of the gauge transformations is obvious: the compatibility condition © is the 
commutativity condition of two linear differential operators 



Qog - g x = , P g - gt = 



L = d x -Q, M = d t -P. 



If the operators L and M commute ([L,M] = 0), then the transformed operators 



L = g- 1 Lg, M = g~ Mg 



(9) 



also commute. Transformations of the form 



L = -h~ 1 L A h, M = -h' 1 M A h 



(10) 
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where L A = —d x + Q tr , M A = —d x + P tr stand for formally adjoint operators and "tr" stands 
for matrix transposition, are also gauge transformations. Gauge transformations form a group, 
the gauge group. 

In the gauge transformations ©, (|1U|) the invertible matrices g, h may depend on x, t, I or 
even be differential or pseudo-differential matrix operators. The Miura and Backhand transfor- 
mations can be viewed as special kind of gauge transformations. 

2.2 Lax representations and Lie algebras 

Equations ©-(JHJ) are naturally defined on any Lie algebra q. Indeed, if Q n ,Pn 6 0, then all 
commutators and derivatives also belong to g. Let q be a finite dimensional Lie algebra over C. 
It is easy to see that only simple algebras are related to coupled nonlinear equations. Indeed, 
according to the Levi-Maltsev theorem ^2]) any finite dimensional Lie algebra q over C can be 
decomposed as a direct sum 

g = K®S 

where 1Z is a solvable radical of g and S is a semi-simple subalgebra. The semi-simple subalgebra 
S, if nontrivial, is a sum of simple subalgebras 

S = es k 

and the following commutation relations hold: 

[n, n) c n , [s k ,n]cn, [s k s n ] = 5 Kn s k . 

If we decompose Q and P in (J2J according to the Levi-Maltsev decomposition theorem 

Q = Qn + Y.Qs k , P = Pn + Y. p s k 

we obviously receive 

d t Q Sk -d x P Sk + [Qs k ,Ps k ] = (11) 
d t Qn - d x P n + [Q n , P n ] + [Qn, £ P Sk ] + £ Qs k ,Pn} = 0- (12) 

We see that equations (jll[) . corresponding to the semi-simple part, can be treated separately 
from the rest of the system and there is no coupling between equations corresponding to different 
simple subalgebras. It is easy to show that, in the radical 1Z, one can choose a basis in which the 
system of equations (|12[1 become "triangular" , i.e. it is an ordered chain of equations such that 
each equation in this chain is a linear equation with respect to its own dependent variable and 
with coefficients and an inhomogeneous part depending on variables of the preceding equations 
and variables corresponding to the semi-simple part. Therefore, the problem of integration of 
the whole system is reduced to the integration of a closed system of the nonlinear equations (|11|) 
and then to the integration of a chain of linear equations 1)12(1 with variable coefficients. 
In this paper we focus our attention on the study of nonlinear equations and therefore we shall 
always assume that the underline algebra q is simple (or semi-simple). 
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2.3 General rational Lax operators 



In this section we define a quite general family of Lax operators with rational dependence on the 
spectral parameter A. Let q be a finite dimensional simple Lie algebra over C and let {qi , . . . , gjv} 
be a basis of q 



N 



[Si,5j)=J2 C ijQr, CT-eC (13) 

r=l 

Let X = {37, CC2, • • •} be a set of independent variables (in the previous example Q we 
used x and i variables). With every variable x k we associate a divisor of poles T k = mi ■ 71 + 
• • • + m n ■ 7 n , i.e. a finite set of points = {71, . . . ,7 n } on the Riemann sphere together with 
their multiplicities mi, . . . , m n £ N. Let C(T k ) denote a linear space of rational functions of the 
spectral parameter A such that the only singularities of these functions are poles at the points 
of the divisor with multiplicity equal or less than the multiplicity of the point. The dimension 
of this linear space dim C^T^) is M k = mi + • • • + m n + 1. Let e\ (I), ... , e k M (I) be a basis in 
C(T k ). 

With x k , Tfc and g we associate the general Lax operator 

d N Mk 

£* = -r-EE^ e S( i )^. (14) 

aXk r=l /3=1 

where u k g are smooth functions of variables X in a certain open domain, and if the set X is 
infinite, we assume that every function depends on a finite number of variables only. A set of 
such functions is denoted by T . 

The Lax operator L k is parametrised by N x M k functions u^g. Any two Lax operators 
L k ,L s form a Lax pair. The commutativity condition [L k ,L s ] = is equivalent to an integrable 
(by the spectral transform method) system of partial differential equations for the functions 



r/3' < V/3' 

In order to write this system of PDEs in explicit form we introduce a basis h^ s (l) , a 
('•'•" ,-' n ''" r ' r - 1 '' ' ; " ] -'"' 



1, ... , Mfc + M s — 1 in /C(Tfc + r s ) and expand elements eg(l),e s gl (I) G -C(r^ + T s ) and the 
products e h Jl)e s g ,(l) € C(T k + T s ) in this new basis 



4( l H' W = Eaii"^" 1 Hf p ?hl s {l) , F*« F|, a , fl*£? e C . 

Thus we have 



(15) 



M fe k Ms a s N M k ,M S 

E^u^-E^af+E E c S .^4^, = o, (i6) 

(3=1 s (3>=1 k i,j=lf3=l,/3'=l 

where r = 1, . . . , N and a = 1, . . . , M k + M s - 1. 

Equation (|16|) is a system of nonlinear PDEs with constant coefficients. The system contains 
partial derivatives in variables x k and x s only, and if we assume that functions u k g , u s T ,g, depend 
on other variables from the set X we can treat these variables as parameters. The system 
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obtained is canonical, it is uniquely denned by the choice of a simple Lie algebra q and two 
divisors and T s (in some cases the divisors may coincide, cf. the Lax pair for the N- 
wave equations pQ). Systems corresponding to a different choice of the basis of g, C(T S ) 
and £(Tfc + r s ) can be obtained from (|16|) by a linear invertible transformations of dependent 
variables. 

The system of equations ()16|) is under-determined, it contains A^x (M^+Mg — 1) equations on 
N x (Mfc + M S ) dependent variables. Indeed, on the Riemann sphere, we have dim£(r/ c + = 
dim£(rfc)+dim£(r s ) — 1 (if we consider a similar setup on an algebraic curve with nonzero genus, 
then the genus would contribute in the accounting for the dimensions of the corresponding linear 
spaces of rational functions on the curve ^Ej)- The difference between the number of dependent 
variables and the number of equations is equal to the dimension of the Lie algebra g. Fixing the 
gauge we finally obtain a well posed problem. 

If we have a family of Lax operators L%, . . . , L\. corresponding to divisors Y\, . . . , Tk, then 
the conditions [Li, Lj] = , 1 < i < j < k yield an over-determined system of nonlinear partial 
differential equations. A general solution and explicit exact partial solutions of this system can 
be constructed by the inverse transform method. 

2.4 The problem of reduction and the reduction group 

The integrable system ()16|) obtained in the previous section is very general. It looks too big 
and may not resemble useful equations but, up to the best of our knowledge, all known 1 + 1 
dimensional systems of integrable partial differential equations are subsystems (reductions) of 
(|lfi|h The problem to find and classify all subsystems of a general system is known as the 
reduction problem. 

Restrictions to subalgebras gCg are obvious reductions of (|16|) . For instance, in the system 
of equations one can consider all matrices to be skew-symmetric and this is compatible 

with the dynamics since the commutator of skew-symmetric matrices is skew-symmetric. With 
this reduction we associate the automorphism eft : a —* —a tr of the Lie algebra sl(N, C). The set 
of all elements of sl(N, C) which are invariant with respect to the automorphisms (j) is obviously 
a subalgebra of skew-symmetric matrices so(N) = {a G sl(N,C) \ a = (p(a)}. A complete 
description of automorphisms of finite dimensional semi-simple Lie algebras is well known (see for 
example |15|). In the case of sl(N, C) automorphisms can be described as follows: for sl{2, C) all 
automorphisms are inner, i.e. they can be represented as (j) : a — > G _1 aG, where G G SL(2,C). 
For sl(N, C) , N > 2 the group of all automorphisms consists of inner <p : a — > G~ 1 aG and outer 
automorphisms V : a -» -H~ l a tr H (G, H G SL(N, C)). 

Let 7£(r) be a ring of rational functions of variable I with poles at points of T and regu- 
lar elsewhere. Fractional- linear transformations of the complex plane /, which map the divisor 
r into itself, induce automorphisms of the ring 7Z(T). We denote by Aut7£(r) the group of 
automorphisms of 7Z(T). Any subgroup of Gr C Autg x Aut7£(r) is a group of automor- 
phisms of q(T) = q<^) c TZ(T). The group Gr is called the reduction group. Elements of the 
reduction group can be viewed as twisted automorphisms, i.e. Lie algebra automorphisms and 
simultaneous fractional-linear transformations of the spectral parameter /. For example, in the 
case sl(N,C), the action of elements of the reduction group Gr can be represented either as 
a(l) -» G- l a{a G {l))G or as a{l) -> -H- l a tr {a H {l))H , where G,H G SL{N,C) and a G (l), a H {l) 
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are the corresponding fractional-linear transformations of the I plane. 

The set QGr (r) = {a £ g(T) \ a = (p(a), \/(f> £ Gr} is a subalgebra of q(T), which we shall 
call Gij-automorphic subalgebra. The restriction of the general Lax operator to the subalgebra 
qg r (T) is a reduction, with reduction group Gr. Restrictions to invariant subalgebras are 
equivalent to require some symmetry conditions for the Lax operators. For example, in the case 
sl(N,C) they lead to symmetry conditions of the form 

L k (l) = G- l L k {a G {l))G , L k (l) = -H- l Li{a H (l))H . (17) 

Comparing with Q and 1)10(1 . we see that the operators are invariant with respect to simul- 
taneous gauge transformations 1 and fractional-linear transformations of the spectral parameter 
I. The reduction group was originally introduced in [H], [H], our recent study of its algebraic 
structure and the corresponding automorphic Lie algebras will be published soon |17j . 

3 ©Tv-Reductions of Lax operators with simple poles, q = s/(2, C). 

In this section we study dihedral reductions of integrable equations corresponding to rational 
Lax pairs with simple poles in the spectral parameter I and Q = sl(2,C). 

The dihedral group Bat is the group of rotations and reflections of the plane which preserve 
a regular polygon with N vertices. It contains iV rotations, which form a normal subgroup 
isomorphic to Zjv and ./V reflections; its order is 2N. If we denote by s the rotation through 
an angle 2ir/N and if r is any of the reflections, then each element of the group can be written 
uniquely either in the form s k or rs k , < k < N — 1. In a more abstract way, the group Djy 
can be defined as the group generated by two elements, s and r, satisfying the identities 

s N = r 2 =id, rsr = s~ 1 . (18) 

In the case N = 2 the group D2 is abelian and isomorphic to Z2 x Z2. 

On the complex plane of spectral parameter A the group ©at can be generated by two 
fractional- linear transformations 

a s : A — > uj\ , o> : A — > A -1 , uj = exp (2m /N) . (19) 

If a divisor V is invariant with respect to transformations (|19|) . it is a union of a finite number 
of orbits of the group D^. The orbit of the group D^v of a point 7 is defined as the set of points 

Dn(i) = R( 7 ), a r <( 7 ) I n = 1, 2, . . . N} . 

When the point 7 is generic, i.e. 7 is not a fixed point of any (nontrivial) subgroup, the 
corresponding generic orbit has 2N points 

Bat(7) ={7, w 7) wV.-.w^-S, 7- 1 , wW-.f^V}- (20) 

As we have already mentioned at the end of section 12.11 there are further generalisations when G and H 
are co-ordinate and / dependent or even (pseudo) differential operators. We do not consider such generalised 
reductions in this paper. 
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Fixed points of the group transformations ()19|) belong to degenerated orbits. There exists a 
degenerated orbit with two elements 

Bjv(O) = {0,oo}, (21) 
and two degenerated orbits with N elements. For odd ./V they are 

Djv(1) = {oj n \n = 1,2,..., AT}, B N (-l) = {-uj n \ n = 1, 2, . . . , N}. (22) 
For even A?" Bjv(1) = Bjv(— 1) and the second orbit with N points is 

B N {u 1/2 ) = {oj n+1/2 | n = 1, 2, . . . , N} . (23) 
In the rest of this section we assume that N is odd. 

With every orbit we associate a linear space of rational functions with simple poles at the 
points of the orbit. Let 7 be a generic point and On (7) the corresponding generic orbit. A 
natural way to construct a basis of rational functions in £(Ojv(7)) is to start with eo = 1 , e\ = 
1/(A — 7) and apply all fractional-linear transformations of the group. As a result we receive a 
basis of £(Bjv(7)) 

eo(Z) = l, e k (l) = , , e N+k (l) = K—r- , k = 1,2, . . . , N. 

to A 7 A — uj K -"-7 

This basis corresponds to the regular representation of the group Bjy. A regular representation 
can be decomposed into a direct sum of irreducible representations. Such a decomposition 
suggests another natural and useful basis of £(Bjy(7)) 

F7 _ . p7 _ A^ + 7^ , 1 + 7^ A^ + 7^ 1 + 7^ 

' 1 - 7^ 1 - 7 Ar A Ar ' 2 A w - 7^ 1 - 7 Ar A Ar ' 

F? - E 1 - F? - E 1 — 



1. 



4fc-l \N _ ryN > 4fc ^ _ jN^N ' ^4fc+l ^ _ \N ' -^4fc+2 \N _ ' 
N-l 



2 ' 

The general Lax operator corresponding to a generic orbit Z?tv(7) is 

2N+1 



L 7 (Z) = ±- - E](l)U] (24) 
3=0 



where Uj are 2x2 traceless matrix functions of independent variables (such as x 7 ). Let L^{1) be 
another Lax operator, corresponding to another generic orbit D^{y) and parametrised by matrix 
functions f/j\ The condition []yy(l), L^il)] = leads to a nonlinear system of 4A?" + 1 matrix 
partial differential equations. We will reduce this huge system to a simple, well determined 
system of six scalar partial differential equations. 

Let a(l) G sl(2, C) (g) c £(B N (~f)), i.e. 

2iV+l 



a(l) = Y E]{l)aj , a,j G *Z(2, 
i=o 
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and consider the following linear transformations 

g s : a(l) -» 5(/i)" 1 a(a s (/))5(/i) , g r : a(l) -> iT^MO)^ . (25) 

where /i is a fixed integer 1 < /i < ^fci that enumerates all different irreducible representations 
and 

s = m = (\" JU), J). <2 6 ) 

The transformations g s and g r generate a linear representation of Bat in the linear space 
sl(2, C) (££) c £(Bjv(7))- It is easy to check that ()25|) generate a subgroup of the group of 
automorphisms of Lie algebra sl(2, C) (^) 7£(Bjv(7))- Let us consider invariant elements of 
sl(2, C) (££) c £(Bjv(t))) i-e- elements a(/) such that 

5(/ l )- 1 a(a,(/))5(/ l ) = a(Z) , Fr l a{a r {l))R = a(l) . 

For every given h the subspace of invariant elements is three dimensional; a basis of this space 
can be written as 



^ " ^ 2 I -1 j ' ^ 2 ~ 1 j ' £[ 



where we fix /i = (JV — 1) /2 since the final result does not depend on the choice of the represen- 
tation. Thus the reduced (invariant with respect to Btv reduction group (|25|l) Lax operator can 
be written as 

L 7 (i) = ^- - - - qz8; , (27) 
where qi,q 2 ,q 3 are scalar functions of independent variables. If 

MO = /" " Pl% ~ P2 £ 2 ~ P3^ , (28) 

is the second operator in the Lax pair, corresponding to another generic orbit Dn(/J>), then the 
condition [L 7 , L^] leads to a system of integrable equations 

qit + a(<M53 - P2<?3)/2 + &(<M?2 - <?3?»3)/2 = 
<?2i - 2/3pig 2 + 47^(^2 + bp 3 )q x = 
q 3t + 2f5p 1 q 3 - (ap 3 + bp 2 )qi = 
Pix + a(g2P3 - P2^)/2 + &(<?2f2 - ^3P3)/2 = 
£>2x + 2ogip 2 - 4:fi N (aq 2 - bq 3 )p 1 = 
P3x ~ 2aqip 3 + ^ N (aq 3 - bq 2 )p x = 

where x = x 7 , t = x M and a, 6, a, (3 are complex constants 

a = (7* - ^Y 1 , b=(l- 7 V)- 1 , 

" - yv _ ^TV + l _ lNfl N ' P- 1 N_ fl N l _ lN/I N ■ 

Clearly equations (|29|) can be rescaled and written in a different form, however they depend on 
parameters 7^ and fj, N , that cannot be removed by simple rescaling. What is surprising is that 
the resulting system does not depend on the choice of the representation. It does not depend on 
N either! Indeed, equations corresponding to N, 7, \x coincide with equations corresponding to 



(29) 
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iVi,7i>Mi> provided the conditions 7^ = 7^ and fi N = fi^ 1 are satisfied. For the case of even 
iV we would receive equivalent equations, reflecting the fact that the corresponding automorphic 
Lie algebras are isomorphic |17j . 

For the orbits Bjv(±1) the bases of £(Bjv(±1)) i n which the representation are decomposed 
in a direct sum of irreducible ones can be written as 

^(0 = 1,^(0 = ^^,^(0 = ^^,^(0 = 3^, fe = i,...,^-. 

The corresponding B^r invariant Lax operators are of the form 

r m d 1 ±F±m ( 1 °\ ±f ^(0 
L±{1) = d^-2 U i F i ®{o -l)- U *[ Ft® 

where again we assume h = (N — 1) /2. The condition [L+, L_] = leads in this case to a rather 
simple system of equations 

U tx- = U 1 U 2 > «aB+ = , = U 2 U 2 > = *4 U 2 • ( 30 ) 

It follows from (J3DJ) that 

«2 = U t U fx„ , U 2 U 2x+ = U l U lx+ , 

therefore we can partially integrate the equations and introduce the variables (j) and 6 

= f(x+) coshcf), U2 = f(x+) sinh^ , u\ = g(x-) cosh 6 , = g{xJ) sinh6> , 

where /(ar+) and g(x-) are two arbitrary functions (which can be set to be equal to 1 by a 
conformal transformation of the independent variables). In terms of </> and equations (j3()[) are 

<fix_ = g(x-)sinh6, X+ = f(x+) sinh0 , (31) 

they are nothing but a well known form of Backlund transformation for the sinh-Gordon equa- 
tion. Indeed, equations (|3*Tj) imply that 

(0 + 4>) x+ x- = g(x-)f(x+) sinh(6» + (j)) , . 

(0 - <t>) x+x _ = g(x-)f(x + ) sinh(# - 0) . ^ > 

Let us take L 7 and L + as a Lax pair. The condition [L 7 ,L + ] = is equivalent to the 
following system of equations 

(g 2 + 53)^2 ( 33 ) 



7 AT_1 

2 



iv_ffe+93)«r ( 34 ) 

9i x+ = 2 (yv_ fj fe + «3)«a ( 35 ) 

47^ + , 2(7^ + 1) + ,o fi , 

q2 *+ = -yv3T^ lU 2 + yv_x W ( 36 ) 

47^ + 2(7^ + 1) + ,„ 7 v 

«3- + = -yv3I^ 2 " yv_i isul . (37) 
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It follows from (O, (JSJ) that 2(u+) 2 + (<u+) 2 = (/(x+)) 2 , where /(x+) is an arbitrary function 
and therefore 

il\ = f(x + ) cos 9 , U2 = V2f(x + ) sin 9 . 
Now equations yield to 



92 + 93 



1-7* 



V2 X7 ' 
In variables 

a o 92-93 



^2(7* + 1) 
equations (|33]) - (|3?|) take the form 



u x+ = f{x + )(cos9) Xl 

v x+ = f(x + )(sm9) Xl (38) 
Ox+x-y = af(x + )usm9 + (3f(x + )vcos9 



where 

( 7 JV_1)2' ^ \~f N ~l, 

Changing the variable x + and rescaling x 7 we can set —* 1, ct a = a/ /3 and /3 — > 1. 

The case of the orbit Bat(O) = {0, co} yields a reducible Lo operator 

d f l\ , 2V- 1 

Lo = d^- p {i-^ o J ' fc = — ■ 

For other values of h the operator is trivial. It is not surprising, since we have imposed the 
condition that the poles are simple. The result become less trivial if we lift this condition (see 

EI)- 

If we take L 7 and Lq as Lax pair, the condition [L 7 ,Lo] = leads to 

1 N 

P x ~,=4qip, 9ix = ^-(93 - 92) , 92xo = -47 qiP, 93x = -49iP- ( 39 ) 

In the new variable u = logp, after simple rescaling and a proper coordinate change, the system 
of equation (|39j) can be written as 



u ^-(u Xl Y 



1ia; 7 xq — 6 

The equation obtained is a well known example of the Liouville type equation ^H]. Considering 
L± and Lq as Lax pair immediately leads to simple C-integrable system of equations. 

For all Lax operators considered in this section the gauge freedom is completely fixed by the 
reduction group. Only scalar I independent gauge transformations commute with all elements 
of the reduction group. 

Equations (|29j) and (|38|) are new integrable systems of equations, to the best of our knowl- 
edge. In particular, (|38|) may have interesting geometrical applications. 
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4 Three examples of D^v reductions, q = sl(N,<C) 



Let us consider a fairly general Lax pair 

L(x,t;X) = d x -X(x,t;X), X = Q + QX + QX' 1 , (40) 

M(x,t;X) =d t -T(x,t;X), T = P + PX + PX' 1 + Q 2 X 2 + Q 2 X~ 2 , (41) 

with Qo, Q, Q, Pq, P, P S sl(N, C). Note that Qo can be always set to zero by a gauge transfor- 
mation. In this setting, the compatibility condition (J2J) yields to the following set of equations 



A 2 : 


Ql- 


--IQ,P) 


A 1 : 


Qt- 


P x + [Q,P ] + [Q,Q 2 ]=0 


A : 


Pq,x 


= [Q,P] + [Q,P] 


-1 . 


Qt- 


P x + [Q,P ] + [Q,Q 2 }=0 


-2 . 


Ql = 


-[Q,P\- 



The system (|42jl is a system of (5(iV 2 — 1)) nonlinear coupled equations for the matrix entries. The 
group of automorphisms of Lie algebra sl(N, C), N > 3 has both inner and outer automorphisms, 
allowing more possibilities for the realisation of the reduction group. Here we consider three 
different reductions of this system with reduction group isomorphic to D^v- 



4.0.1 Case 1: inner and outer automorphisms. 

Let us consider two transformations 

s : L(l) ^ S' 1 L[ul) S , r : L(l) ^ -L A (l/l) , (43) 

where S is a N x N matrix given by Sij = 5ijio N ~' t , with u = e 2t7T ^ N , and where li L A " stands for 
formally adjoint operator L A = —d x + X tr . Observe that, neglecting the spectral parameter I, 
1)43(1 are nothing but two automorphisms of the algebras sl(N,C); in particular, the first one is 
an inner automorphism, while the second is outer ^]. They satisfy (|18|) and therefore generate 
the dihedral group Dat. 

Let us now require that both operators L and M are invariant under (|43jl. That leads to 
algebraic constraints on the matrices X and T 

X(l) = S- l X(ul)S, X(l) = -X tr (l/l), 

T(l) = S~ 1 T{ul)S, T{1) = -T tr (l/l), { > 

where "tr" stands for matrix transposition, which imply 

^o = 0, 

Q = -Q tr , P = -P tr , 

Q = q(x,t)A, P = p(x,t)A, 

where qij = qi(x,t)5ij and = pi(x,t)5ij are diagonal matrices and where A is the shift 
operator Ay = 5ij-i. Here and hereafter all indexes are counted modulo N. As a consequence, 
(|42|) reduces to the following system of 2(iV — 1) nonlinear coupled equations 



Qt-P x + [Q,Q 2 } = o 
Ql = [Q,P] 



(45) 
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or, in components, 

qu ~ Pix + qiq 2 i+ i - q 2 i-iqi = ( 46 ) 
qiPi +1 - Piq i+ i - {qiqi+i) x = . (47) 
Solutions of (|47l) can be parametrised by new variables itj and Vi (similarly to jSJ); indeed, let 



qi = exp(nj) , Pi = Vi exp(uj) , (48) 

then 



A /2mod(i-r- 1,N) + 1-N\ , 



r=l x 

In the new variables (|4fi|) reads 



uu - v ix - Ui x Vi + exp(2-ii i+ i) - exp(2«j_i) = 0, i = l,...,N. (50) 
In the case N = 3 and Q 3 = I (or Y\ qi = 1)> equation (|4"5|) can be solved explicitly 

P = \[Q,QQxQ\. (51) 

In components 

Pi = ^<?l(<?2<?3x - fe?3) , P2 = ^q%(qmx ~ q3xqi) , P3 = ^(qmx ~ <?Lr<72) • 

Hence, substituting into (|4f)]) and rewriting the equations in terms of Ui variables we obtain 
uu = ^{u 3xx - u 2x x) + t^u 1x (u 3x - u 2x ) ~ exp(2u 2 ) + exp(2u 3 ) (52) 
and cyclic permutations of the indexes 1, 2, 3. 



4.0.2 Case 2: inner automorphisms. 

Let us now turn our attention to a different symmetry conditions 

L(l) = S- 1 L(wl) S , L(l) = R- 1 L(l/l)R (53) 

where Rij = 5i,N-j (all indexes are counted modulo N). Conditions (|53jl are twisted inner 
automorphisms. In this case 

(Po)ij = PoA,j , 

Q = RQR, P = RPR, 

Q = q(x,t)A, P = p(x,t)A. 



Hence (|4"2*|) reduces to 



Qt-Px + [Q,P } + [Q,Q 2 } = o (54) 

P 0x = [Q,P] + [Q,P] (55) 
Ql = [Q,P], (56) 
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or, in components 

q%t = Pix + poiQi - QiPoi+i + qiQi+iqN-i-2 - qN~iqi-iqi (57) 
Poix = qiPN-i-i - PiqN-i-i + qN-iPi-i - PN-m-i (58) 

qiPi+i -piqi+i - {qiqi+i) x = 0, i = l,...,N. (59) 
For JV = 3we can use the result (|51|) for P to solve (|55|) and find 

Po(x, t) = q 2 (x, t)q 3 (x, t) diag {1 ,1,-2}. (60) 
In U{ variables (|48j) . we obtain 

Ult = k(u3xx ~ U 2x x) + kuix{u3x ~ U 2x ) 

U2t = 3 {u\ xx ~ u 3xx ) + zU 2 x{uix ~ u 3x ) - exp(2ui) + 4exp(u 3 + u 2 ) (61) 
u 3 t = 3 {u 2x x ~ uixx) + lu-3x{u2x - uix) - 4exp(ii 3 + u 2 ) + exp(2tii) . 



4.0.3 Case 3: symplectic automorphisms. 



In even dimensions N = 2n we can consider a third case given by 

L(X) = S- 1 L(uj\)S , L(A) = -J- 1 L A (1/X)J , 



(62) 



where J is the symplectic matrix J = (V„). From (|62|) it follows that Pq is a diagonal matrix 



while 



Q 



Po = 


diag{p i(^,i), • 




,P0n(x,t) 


-Pai(x,t), 


• ' ' 1 


-POn 


{x,t)}, 


9i 


... 








( 


Pi 











q 2 ■■■ 














P2 








'•. 




P 



















' • • 92n-l 

















' ' • P2n 


q2n o 


... 


) 






\ -Pin 











(63) 



and 



Q = -J- 1 Q tr J p 



-J- l P tr J . 



The equations for this case read 

q%t = Pix + qiPoi - poi+m + q 



n+%~\qi-\qi — qiqi+iqn+i+1 

POix = —qiPn+i + q-n+iPi + qi-lPn+i-1 ~ q-n+i-lPi-1 

qiPi+i - Piq i+ i - (qiq i+ i) x = , i = l,...,2n. 



(64) 
(65) 
(66) 



For n = 2 we have [Q , Q 2 ] = and equation (|64j) simplifies. In terms of Uj variables (|48|) 
and ()49|) the system becomes 

Wu + \{w2xW3x) = 

w 2t + w 3xx + l(w 3x wi x ) - 2l 2 = (67) 
w 3t ~ w 2xx + \ {w 2x wi x ) - 2h = 

hx = 2w 2x exp(-wi) , l 2x = -2w 3x exp(^i) , (68) 
where w\ = u\ + u 3 , w 2 = u\ - u 3 , w 3 = u 2 - u 4 and where h = (p m + p 02 ), l 2 = (p i - p 02 ). 

System ()67(l is new, to the best of our knowledge. Moreover, for TV" general all systems of 
equations obtained in this sections can be regarded as new examples of lattice equations. 
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5 Conclusions and remarks 



All properties of integrable equations are encoded in their Lax representations. Therefore, the 
description of the variety of integrable equations and their reductions can be pursued starting 
from the theory of Lax operators. In particular, the problem of reductions can be studied impos- 
ing symmetry conditions on the Lax operators, reducing the corresponding nonlinear systems of 
equations to smaller subsystems. Symmetries of Lax representations form a group - the reduction 
group |H1, 0, JU], ^7j. The reduction group approach has been used to describe Lie-algebraic 
reductions of the iV-wave equation to find Lax representation for a number of new Non- 
linear Schrodinger type equations ^H! an d to build explicit solutions for the Landau-Lifschits 
equation [20] • In this paper we have revisited the reduction group approach, illustrated it by a 
number of examples and motivated our further study. There are several natural directions for 
development. 

The equations obtained in the previous section are integrable for any N. Assuming N —* oo 
and taking continuous limits one can find 2 + 1 dimensional integrable equations (and corre- 
sponding Lax representations - all structures, such as symmetries, conservation laws, etc. can 
be recomputed in these limits). There are several ways to take a continuous limit, the result 
depend on the balance of non-linearity and dispersion. For example, the Kadomtsev-Petviashvili 
equation 

U T = -Uyyy - D^tllXx) +6UUy . (69) 

can be recovered from Q46j) - Q47JI as a continuous limit if we assume q = exp(/i 2 u), where h is the 
lattice step, and perform a Galilean transformation. On the other side, performing rather naive 
expansions gy-i = q ± hq y + 0(h 2 ) and p%±\ = p ± hp y + 0(h 2 ), we would receive, after proper 
rescaling, the hydrodynamic type equation 

ut = DyiiD^Dx) 2 u + \{D- X D X u) 2 + e < 2u >} . (70) 

Similarly, equations (|57|) -(|59 |) yield to the following systems of equations in 2 + 1 dimensions 

u t = D y {{D y l D x ) 2 u + \(D~ X D X uf + D~ l D y [e^ D^D^u + v)] + e^} 

vt = -DyKD^Dx) 2 v + \{D~ l D x v) 2 + D- l D y [e^ D^D x (u + v)\ + e^} , 1 j 

while and (jBU)-(jSni) yield to 

u t = D y {{D- x D x ) 2 u + \{D y x D x u) 2 + D- l D y [e^ D^D^u - v)} + e^} , . 

v t = D y {(D y 1 D x ) 2 v + \{D y x D x v) 2 + D~ l D y [e^ D^D x {u - v)\ + e^)} . 1 ' 

Notice that ()7Uj) can be obtained from (|72|) by setting v = u. 

A simple modification of the procedure enables us to study Lax operators with non com- 
mutative (matrix) entries. For example, we can treat elements of sl(2M, C), as 2 x 2 matrices 
with M x M matrix entries. The corresponding reduced systems can be viewed as a system of 
equations for non abelian variables. The following system of equations 

qot - Pox + [qo , Po] + [qi , Pi] = 
qit + [qi , po] + {q2 , P2>/4 = o 

q2t + [qa , Po] + {qi , P2} = (73) 

Pla: + [Pi , qo] + {q2 , P2}/4 = 

P2x + [P2 , qo] + {q2 , pi} = o 
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where the variables q, , Pi are elements of a non commutative free algebra (or matrices of any size) 
and {qj , pj} = qiPj+p^qi, is a nonabelian generalisation of Q31|) (the Backlund transformations 
for the sinh-Gordon equation). In the non abelian case we have to fix the gauge freedom further 
in order to make equations (|73|) well determined. 

A very important issue is the solution of nonlinear integrable models. The reduction group 
proves to be not only a very useful tool to find new integrable equations and classify Lax pairs 
but also a necessary instrument in this context. Indeed, without the constraints imposed by the 
reduction group on the spectral data it is not even possible to formulate the inverse problem, 
which would lead to explicit solutions of the integrable equation. The general set up of the 
correspondence between reduction groups and analyticity properties of the spectral data is one 
of the next issues on our research agenda. 

Another challenging problem is the study of automorphic Lie algebras in a pure algebraic way. 
They can be always linked back to Lax operators, Baxter's i?-matrix equations, etc. Moreover, 
the problem of a complete description of rational automorphic Lie algebras seems to be feasible. 
Our optimism is based upon a simple group-theoretical observation and a remarkable the- 
orem of Felix Klein the complete list of finite groups of fractional-linear transformations is 
given by the cyclic group "Lm, the dihedral group Bat and the groups of symmetry of Plato solids, 
i.e. the tetrahedral group, the octahedral group and the icosahedral group. It would be interesting 
to generalise the theory of automorphic Lie algebras to the cases of elliptic and higher genus 
algebraic curves. We believe that automorphic Lie algebras will find applications far beyond the 
theory of Lax operators. 
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